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Abstract: We study the evolution of a two component fluid consisting of "blue" 
and "red" particles which interact via strong short range (hard core) and weak long 
range pair potentials. At low temperatures the equilibrium state of the system is one 
in which there are two coexisting phases. Under suitable choices of space-time scal- 
ings and system parameters we first obtain (formally) a mesoscopic kinetic Vlasov- 
Boltzmann equation for the one particle position and velocity distribution functions, 
appropriate for a description of the phase segregation kinetics in this system. Further 
scalings then yield Vlasov-Euler and incompressible Vlasov-Navier- Stokes equations. 
We also obtain, via the usual truncation of the Chapman-Enskog expansion, com- 
pressible Vlasov-Navier-Stokes equations. 



1. Introduction. 

The process of phase segregation in which a system evolves from an initial unstable 
homogeneous state into a final equilibrium state consisting of two coexisting phases is 
of continuing theoretical and practical interest [GSS], [FLP], [L]. Such a process occurs 
whenever the system, which is initially at values of the thermodynamic parameters, say 
temperature To and pressure po, corresponding to a single homogeneous phase has its 
parameters "suddenly" changed to new values, say T and p, at which there is a coexistence 
of phases. 
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This happens, for example, when an alloy is 'quenched' from a high temperature melt 
or solid to a low temperature solid state by sudden cooling [GSS]. After such a quench the 
system finds itself in an unstable (or metastable) situation, as far as the spatial concen- 
trations, which have not been able to adjust rapidly enough to the "sudden" quench, are 
concerned and domains of the two equilibrium phases form and start growing in time. This 
proceeds until there are "in the final state" only two regions of pure equilibrium phases 
separated by an interface. Since the kinetics of the domain growth have a profound influ- 
ence on the properties of the alloy, this problem has been and continues to be extensively 
studied both theoretically and experimentally [GSS]. For such alloy systems the segrega- 
tion process takes place mainly through the (anti) diffusion of the two components — from a 
uniformly mixed state to a demixed one. There are no macroscopic matter or energy flows 
since the system is a solid and has a high heat conductivity which keeps the temperature 
equal to some constant ambient value. The only relevant conserved quantities are therefore 
the particle numbers of the two components and the macroscopic equations describing the 
process are fairly well established: these are the well known Cahn-Hilliard equations [CH] 
and variations on them. We refer the reader to reviews on this subject [GSS]. 

The situation is much less clear for phase segregation in fluids where macroscopic flows of 
matter and heat are important. There are now additional conservation laws for momentum 
and energy and there is no general consensus even on what hydrodynamical equations 
are most appropriate for describing the macroscopic evolution of the system [S], [OP], 
[AG A] . In particular, it is not clear which is the correct coupling between the Cahn-Hilliard 
equation for the order parameter and the Navier-Stokes equation for the fluid velocity. 

To make a start on the mathematical analysis of such processes we investigate a model 
binary fluid introduced in [BL] where the process of phase segregation was studied numer- 
ically. 

In the present work we derive general equations appropriate both in the one phase and 
in the coexistence region. In part II we consider applications to the segregation process 
including an analysis of new numerical results. Many of our discussions here will be semi- 
heuristic. In particular, we will not go into detail about the domain of validity of the 
technical conditions necessary for the rigorous mathematical establishment of the results. 
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The model we study is composed of two types of particles, call them red and blue. 
There are N r red and Nb blue particles in a cubic box of volume A = L d ; we will generally 
consider d = 3 and use periodic boundary conditions. The particles all have unit mass and 
hard core diameter a. Particles of different kind also interact with each other through a 
long range pair potential of the Kac type, having a range i and a strength Ag. By properly 
choosing At, we obtain, in the limit I — > oo, a system whose equilibrium properties are 
described by a mean-field type phase diagram exhibiting a demixing phase transition for 
temperature T < T c [LP]. 

This transition is essentially independent of the hard core size a and the dimension- 
less microscopic particle densities p r a 3 and p&a 3 can therefore be arbitrarily small in the 
demixed phases. This means that we can have a situation in which, at least in principle, 
the whole phase transition is well described by a Vlasov-Boltzmann type of kinetic equa- 
tion. We will in fact see that we can, by suitably scaling space and time and the densities, 
obtain, at least on the formal level, a set of nonlinear Vlasov-Boltzmann (VB) equations, 
describing the evolution of the one particle distribution functions f a (q, v,t), a = r, b. 

The VB equations we derive are of a form similar to ones conjectured for a one com- 
ponent fluid with hard cores and an attractive long range interaction [DS], [G]. Such a 
system however requires the hard cores for stabilization against collapse [LP] and pa 3 is 
greater than 1/3 in the liquid phase. It is therefore not clear that a VB equation is an 
appropriate kinetic description of such a liquid-vapor transition. This is the motivation 
for introducing the binary model we consider here. 

We discuss the scalings necessary to go from a microscopic Hamiltonian description of 
the time evolution to the VB equations in Section 2 leaving a formal derivation, in the 
spirit of Cercignani [C] and Lanford [Lan], to Appendix C. The equations themselves are 
of the same form as those used in [BL] for the kinetics in the coexistence regime of this 
system. Their numerical results for the time evolution and the analysis of the stationary 
states showed that these VB equations for the one particle distributions f a indeed lead to 
the phase segregated state expected from purely equilibrium considerations. 

While the mesoscopic description in terms of the one-particle distribution functions 
is a great simplification compared to the full microscopic representation, it is still more 
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complicated than the macroscopic theory that treats the binary system as a continuum 
with well defined local density p(x,t), concentration difference cp(x,t), velocity u(x,t) 
and temperature T(x,t). The derivation of hydrodynamic equations from the Boltzmann 
equation (which one expects to be structurally of the same form as those describing dense 
binary fluids) is closely related to the problem of finding approximate solutions of the 
Boltzmann equation. The reason for this is that the fluid dynamic variables are defined 
and change on space and time scales which are very large when measured in units of the 
mean free path and mean free time between collisions, i.e. the kinetic or mesoscopic scale. 
Therefore, it can be expected that the system will reach a state close to local equilibrium 
in a macroscopically very small time interval, meaning that f a {x,v,t) should stay close 
to local Maxwellians, with parameters p a , u and T, which change slowly on the kinetic 
scale. The big disparity between the kinetic and hydrodynamic scales suggests looking for 
a solution of the Boltzmann equation as a series expansion in the scale parameter which 
is the ratio of these two scales. Many rigorous results in this direction have been obtained 
in recent years, especially for the Euler (E) and the incompressible Navier-Stokes (INS) 
equations. The situation is less satisfactory in the case of the compressible Navier-Stokes 
(NS). This is a consequence of the fact that while the E and INS equations correspond 
to well defined scaling limits, in which the mean free path goes to zero, there is no such 
scaling limit for the NS equations as can be seen from the fact that these equations are 
not invariant under scaling [DEL]. 

Having obtained the VB equations we turn to the derivation of hydrodynamic equations. 
The results available for these equations are fewer than for the Boltzmann equation. In 
Sect. 3 we present a rigorous derivation of the Vlasov-Euler (VE) equations for this system, 
which differs from the usual Euler equations by the presence of self-consistent forces coming 
from the Vlasov terms. We do this by adapting to this case the method of Caflisch [Ca80] , 
i.e. we prove that the Hilbert expansion is asymptotic, by showing that the remainder at 
any order is finite in a suitable Sobolev norm. 

We then consider in Sect. 4 and 5 a modified Chapman-Enskog expansion of the kind 
considered by Caflisch [Ca87] and show also in this case that the remainder at any order is 
finite in the same Sobolev norm. The term of zero order in this expansion is a Maxwellian 
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with parameters solving a set of dissipative new PDE's, the Vlasov-Navier-Stokes (VNS) 
equations, where, beyond the usual terms present in the compressible Navier-Stokes equa- 
tions, there are diffusive terms coming from the presence of the self-consistent force. In 
particular, the equation for the concentration can be put in the form of a gradient flux 
of an energy functional [BELMII] which is similar to an exact evolution equation derived 
for a microscopic model of a binary alloy. The latter has been proven to yield the same 
late time phase segregation behavior as the Cahn-Hilliard equation, [GL96] , [GL97] . Both 
Vlasov-Euler and Vlasov-Navier-Stokes have non trivial stationary solutions with the same 
solitonic profile as in the BV equation. 

Finally in Sect. 6 we consider the incompressible regime for these equations and derive, 
under suitable initial conditions and scaling, a set of PDE's with dissipative terms involving 
a force linear in the concentration (they are essentially the linearization of the analogous 
terms in the compressible equations around a constant concentration and density profile) . 
Above results all rely on the crucial assumption that the initial value problems for the 
hydro dynamical equations have a unique smooth solution at least on some macroscopic 
time interval. We do not discuss the technical conditions which ensure the existence of 
such solutions. 

2. Vlasov-Boltzmann equation for a binary mixture. 

We consider a system of N r red particles with positions and velocities v\, i = 1, . . . , N r 
and Nj, blue particles with positions and velocities v\, i = 1, . . . , Nb, in a 3-dimensional 
torus A, interacting via two body forces. N = N r + N b is the total number of particles. 
The potential energy is 



nff & r ;a &j = ^ EEE w-^i) 



aj=/3 i=l 3=1 



(2.1) 




a,f3 i=l j=l 
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where a, (3 = r, b, Ue is the long range potential 



Ue(r) = U{ r -) 



(2.2) 



for some bounded, smooth non-negative function U on IR + . The factor Ai is the intensity 
of the long range interaction to be suitably chosen to get a mean field type of behavior 
when I becomes very large compared to the interparticle spacing [LP]. The potential W a 
is the formal hard core potential 



In other words the particles are hard spheres of diameter a interacting by elastic col- 
lisions which are color blind and by a weak repulsive long range force between particles 
of different species. The total number of particles of each species as well as the total 
momentum and energy are invariant during the evolution. 

Choosing the size of A to be I (or some constant multiple thereof) there are two charac- 
teristic length scales for this dynamics: a, the range of the hard core potential and £, the 
range of the Kac potential. We can consider a third length, which depends on the density, 
the mean free path A defined by the relation 



The kinetic limit arises when there is a large separation between a and A, corresponding 
to a low density (N/£ 3 ) situation. To obtain a kinetic limit we send N and £ to oo while 




if r < a 
otherwise. 



A = 



Na 2 ' 



a is fixed, say 1, in such a way that X/£ is finite (N ~ £ 2 ) and assume initial data almost 
constant on regions of size ~ A. We denote by 




1 



A 




The kinetic equations will be obtained in the limit 5^0, assuming 



A e = 7 3 5 2 , 
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meaning that Ag is proportional to 1/N. A further limit 7^0 will provide the hydrody- 
namical limit to be discussed later. 

In kinetic coordinates q, that is qf = 6£f, for a = r,b and i = 1, . . . , N a and kinetic 
time t = 5r m , r m being the microscopic time, the equations of motion for the system are, 
for a = r, 6 and i a = 1, • • • , N a 



dqf 

a = x) 



dr 

r r a *e ( 2 -3) 

^=7^E^(7lC-<l)(l-M 
.7/3=1 



m 

, N „ td3N 



T N = {(q 1 ,v 1 ,...,q N ,v n ) e A w x IT ||ft-^|>5, i ^ j}, 

where K(\x — y\) = —(VU)(\x — y\), N = N r + Nb, and we use the notation qi,Vi, 
i = 1, . . . , N when the color is irrelevant. When two particles are in contact (namely at 
distance 5) they undergo an elastic collision regardless of their color. We neglect the event 
that more that two particles are in contact because it has vanishing Lebesgue measure. 
Hence the evolution (2.3) is defined only almost everywhere. 

Note that when N ~ £ 2 then the mean force on each particle on the kinetic scale (2.3) 
is of order unity. This is the reason why our original choice of the strength Ag of the 
potential in (2.2) was like £~ 2 rather than £~ 3 as in the usual case [LP]. 

With this scaling we can get, at least formally, in the limit 5 — * the Vlasov-Boltzmann 
equation for a binary mixture of hard core particles interacting via a weak long range 
potential. A formal proof of this is given in Appendix C. The rigorous proof would require 
the extension of the Lanford argument to this case, an extension that is not obvious because 
the Vlasov part is not well controlled in the Lanford norms. 

Even if we have discussed the derivation of the Vlasov-Boltzmann equation only for hard 
spheres, from now on we consider the Vlasov-Boltzmann equations in full generality. The 
function f r (q,v,r) (resp. f b (q, v,r)) is proportional to the probability density of finding 
a red (resp. blue) particle at q G O C K 3 , with velocity v G K 3 at time r > 0. We notice 
that the relation between the / Q 's and the microscopic one particle densities Pi(£,v,T m ) 



(normalized to iV a ) is given by 



r(q,v,T) = \im6- 1 p?(6- 1 q,v,6- 1 T). 

o — >0 

The functions f r and f b are positive and normalized to 7 -3 for any value of r. They 
are solutions to the equations 

o T f r + v ■ v q r + F r ■ v v r = j(r , n + Ar, A 

(2.4) 

d T f b + v ■ VJ b + F b ■ V v f b = J(/ 6 , f b ) + J(f b , f). 

The Vlasov force acting on each particle is of the Kac type, meaning that for any 7 > 0, 
the forces are conservative non local forces with range 7 _1 defined by the position 

F a (q,r) = -V q [ dq'^U^q-q'Dn^q^r), a = r,b, a ^ (3 (2.5) 

with {7(|<z|) a smooth, non negative function of compact support and n r , n b are the rescaled 
spatial densities of the red and blue particles: 

n a {q,T)= [ dvf {a) {q,v,r), [ dq n (a) (q, r) = 7 " 3 . (2.6) 
Jtr 3 Jo. 

For any positive functions / and g, J(f,g) denotes the effect of the collisions of particles 
distributed according to g on the distribution /. Its expression is given by 

J(f,g)= [ dvj dub(\v-v^u)[f(v')g(v:)-f(v)g(v,)]. (2.7) 

JjR 3 JS 2 

Here 6(]v|,o;) is the differential cross section of the short range interaction, u G S% is the 
impact parameter and v', v'* are the incoming velocities corresponding to an elastic collision 
with outgoing velocities v,v* and impact parameter u. We assume the Grad's (see [Gra]) 
angular cutoff condition that 6(|v|,a;) is a smooth function growing at most linearly for 
large \v\, i.e. 6(|u|,a;) = \v\ a h(u) with < a < 1 and h a smooth bounded function on S2. 

An important property of the collisions is the entropy production inequality: let 

M a = [ dvJ(rj a )\ogf a , a = 1,2, 
Jr 3 



M a ,p = [ dvJ(r,f)logr, a, 0=1,2. 

Then A/" Q 's as well as A/1,2 + A/2,1 are non negative. Moreover iV a vanishes as usual if and 
only if the / a, s are Maxwellians: 

f a = M(n a ,u a ,T a ;v), a = 1,2 

with 

M(n, u, T; v ) := e"^ ^ ^ . (2.8) 

Furthermore A/1,2 + A/2,1 vanishes if and only if the two Maxwellians have the same local 
temperature and mean velocities: 

u a = u, T a = T, a = 1, 2. 

This implies that the only solutions of the equations 

J(/i,/i) + J(/i,/ 2 ) = 0, 

J(/2,/ 2 ) + J(/2,/l) = 

are Maxwellians with the same mean velocity and temperature. General arguments suggest 
that all the stationary solutions of equations (2.4) will be Maxwellians with u = 0, T(q)=T, 
and densities satisfying the equations 



Tlogn 1 (q)+ J dq' 1 3 U( 1 \q-q'\)n\q , ) = C 1 
T\ogn\q) + J dq'^U^q - q'\)n\q') = C 2 



(2.9) 



Beyond the spatially constant equilibria, there may be other spatially non homogeneous 
solutions. For example, by prescribing the boundary conditions in one dimension 

lim n a {z) = n+, 

one gets at small values of T a solitonic solution describing the interface profile [BL]. We 
shall leave a discussion of this part for [BELMII] and focus here on deriving macroscopic 
equations for the evolution of the conserved quantities. 

9 



Before closing this section, let us define 

f(q,v,r) = ^[r(q,v,T) + f b (q,v,r)] 

as the density of finding a particle at q with velocity v at time r, independently of its color. 
Moreover, we set 

<j ) (q,v,T) = ^[r(q,v,r)-f b (q,v,r)}. 
The system (2.4) can be written in the following equivalent form: 

drf + v- V q f + 2F ■ V v f + 2W • V„0 = 4J(/, /), 

(2.10) 

d T (f> + v ■ V q( f> + 2F ■ V V <P + 2W ■ V v f = 4J(0, /), 

where F = F r + F b : W = F r — F b . We can absorb the numerical factors by redefining U 
as U/2 and b as 6/4 so obtaining 

d T f + v ■ V q f + F ■ V v f + W ■ V„0 = J(f, /), 

(2.11) 

d T <j) + vV q <P + F-V v <P + W ■ V v f = J(0, /). 



3. Compressible Hydrodynamics. 

We are interested in the behavior of the system on the macroscopic scale. To this end we 
introduce a scaling parameter e representing the ratio between the kinetic and macroscopic 
space units and, for any t > and x G eVL we set 

r = e~ 1 t, q = e~ 1 x. 

We assume that at time zero the densities vary slowly on the microscopic scale f l {q, v, 0) = 
f(eq,v,0) and look for solutions of (2.11) such that 

f(q,v,r) = f(eq,v,er), i = r,b 

with f l smooth functions on eQ x IR 3 x IR + . For the force we have 

F r (e- 1 x,£-H) = - £ ^xf s dx '{~Y U (^)\x-x'\\n b (x',t), 
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and a similar relation for F b . Therefore, if we assume 7 = e, also the forces are slowly 
varying functions and the / satisfy the following system, where we remove the "tilde's" 
because in the sequel we shall always consider only the macroscopic variables: 



d t f + v-V x f + F-V v f + W- V v <j> = e- 1 J(f, /), 
d t 4> + v ■ V x + F ■ V v (p + W ■ V v f = e- 1 J(0, /). 
We shall use the notation 



(3.1) 



F = K®f, W=-K ©0, (3.2) 

where, 

K(x) = -V x U(\x\) (3.3) 

and, for any function g we set 

(K© g)(x,t) d = [ dx'K(\x-x'\) [ dvg(x',v,t). (3.4) 
Jn Jr 3 

We will show that the solution of the system (3.1) is close for e small to the local 
equilibrium with parameters p^ 1 ) , , u and T satisfying the following set of hydrodynamic 
equations: 

d t p + v • h = 0, 

d t (p + V • {(fu) = eV ■ (DQ), 

(3.5) 

pD t u + VP — pK * p + (fK * if = -eVff, 
3 

-pD t T + PV ■ u = 6V(kVT) - so : V« - eK * <p ■ DQ. 

Here p = p^ + p^ is the total density, if = p^ - p (2) , P = pT 

D t :=d t + u- V, 

a ■= -i/(V« + Vu ] 1 V • u) 

3 

Q:=V^ + -^-(p 2 -<p 2 )K*<p. (3.6) 
P P T 

is the adjoint of the matrix Vw, a : Vu = Tr (crVw), I is the unit matrix, v and 

D/p are the viscosity and the diffusion coefficients and k is the heat conductivity. These 
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are computed from the VBE. The above equations, with e = will be referred to as the 
Vlasov-Euler equations (VE). We assume that the initial value problem for such equations, 
with suitable initial data, has a sufficiently smooth solution at least on a time interval [0, t\. 
Under such conditions we will prove in the next section and in Appendix A that the solution 
to the VBE for the binary fluid, under the Euler scaling, converges to the Maxwellian local 
equilibrium with parameters satisfying the VE equations in the interval [0, t\, with an error 
of order e (Proposition 4.1 and Corollary 4.2). 

When e > 0, the above equations will be referred to as the Vlasov-Navier-Stokes equa- 
tions (VNS). In Section 5, using also the arguments of Appendix A we will show that 
their solutions provide an approximation up to the order e 2 to the solutions of the VBE in 
the Euler scaling, provided that the initial value problem for such equations has suitably 
smooth solutions as before. The precise statement is given in Proposition 5.1 and Corollary 
5.2. 

Like for the usual Navier-Stokes equations, which are frequently and successfully used 
with e = 1 in physical and engineering applications, although their derivation is restricted 
to small values of e, we will consider the VNS equations with e = 1 and analyze some of 
their properties in [BELMII]. In order to get diffusive effects as sharp limits of the VBE, 
it is necessary to go to the parabolic scaling where r = e~ 2 t and consider simultaneously 
a low Mach number situation. This will be discussed in Section 6. 



4. Euler limit. 

We outline the proof of the convergence of the Vlasov-Boltzmann system to the VE 
equations. The proof will be completed in Appendix A. We fix the Maxwellian M(p, u, t) 
with p, -u, T possibly depending on space and time and denote 



Cf = J(M,f) + J(f,M) 
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(4.1) 



and 



r/ = j(/,m). 



(4.2) 



Moreover, we set 



Q{f,9) = l[J(3,f) + Jtf,9)] (4-3) 



It is easy to check that, as for the one-component Boltzmann equation, 

Cf = iff f = M Xa , a = 0,...,4, (4.4) 

where 

Xo = 1, Xi = Vi, i = l,..., 3, X4 = v 2 /2 (4.5). 
Moreover, along the same lines one gets 

17 = iff f = aM, a G M. (4.6) 

We shall try to solve (3.1) following [Ca80], in terms of a truncated Hilbert expansion 
of the form 

K 



K 



(4.7) 

K 



n=0 

with suitably chosen positive integers K and m. The functions f n and 4> n are computed 
using a Hilbert expansion and the remainders Rf and are defined as the difference 
between the solution and the truncated expansion. 

We substitute in (3.1) the formal power series 

oo oo 

/ = Z> n /»> 0=$>>n, (4.8) 

n=0 n=0 

oo oo oo oo 

F = ]T e n F n = ^ K ®fn, W=J2 £UW " = J2 £UK ® tn, (4-9) 

n=0 n=0 n=0 n=0 
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and denote by D t the time derivative along the trajectories: 



D t = d t + v ■ V 3 



We have: 



00 

£ _1 Wo, fo) + $> n [ 2 ^+1) + S n ] = 0, 



n=0 



(4.10) 



00 

e- 1 J(<h, fo) + J2 £U [ J (<^+i, fo) + J(<fo, fn+l) + T n ] = 0, (4.11) 



n=0 



where 



S n = Q(fhJh>)~ V [F h -V v fh>+W h -V v (f> hf ]-D t f n . (4.12) 



(h,h'): h,h'>l 
h+h'=n+l 



(h,h'): h,h'>0 

h + h'=n 



J(<ph, h>) - i Fh • v ^ h ' + Wh • v -^'] - ( 4 - 13 ) 

(h,h'): h,h'>l (h,h'): h,h'>0 

h+h'=n+l h+h'= n 

In order that the formal series solve (3.1) the coefficients have to satisfy the conditions: 



Wo,/o) = 0, 
J(<h,fo) = 



(4.14) 



and, for any n > 0, 



(4.15) 



2Q{f ,f n+ i) + S n = 0, 

J((j>n+l,fo) + J((f>0, fn+l) +T n = 0. 

As remarked in the previous section, the first of the conditions (4.14), implies that fo is a 
Maxwellian with parameters depending on x, t: 



f (x, v, t) = M(p(x, t),u(x, t),T(x, t);v) := M(v). 



(4.16) 



Moreover from the second eqn. of (4.14) we get 

(p(x,t) 



(f) (x,v,t) = ' M(p(x,t),u(x,t),T(x,t);v), 
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for some suitable function ip(x,t). Using (4.1) and (4.2) we can write (4.15) as 

Afn+l = — «5 n , 

(4.17) 

T0n+1 = -J(<j>0, fn+l) - %i- 

Since <S n and T n only depend on the fk and 4>k f° r k < n, we have first to solve the first 
eqn. and then, once / n +i is determined, we solve the second one for (p n +i- 

In order to check the solvability of these equations we introduce the Hilbert space H of 
measurable functions on IR such that the scalar product 

(f,g)= [ dvf(v)g(v)M-\v), (4.18) 

is finite. In this Hilbert space the operators C and V are densely defined and symmet- 
ric. Moreover, the null spaces of C and V are the five-dimensional subspace spanned by 
{Mx Q , a = 0, . . . , 4} introduced in (4.5) and the one-dimensional space spanned by Mxo 
respectively. We denote by V and /C the projectors on such subspaces and by V 1 - = 1 — V 
and JC 1 - = 1 — K the projectors on their orthogonal complements. From the properties of 
C and T it is immediate to check that 

VC = 0, KT = 0. (4.19) 

Both C and V are non positive and there are positive constants 5 and 5' 

(/, Cf) < -5||^/|| 2 , (/, r/) < -5'W^fW 2 . (4.20) 

Moreover 

C=-v + K, 

(4.21) 

r = -u + e, 

where 

v(x,v,t)= dv*b(\v-v*\)M(p(x,t),u(x,t),T(x,t);v*), (4.22) 
is a strictly positive function such that 

fo(l + M) CT < < ^i(l + |v|) CT (4.23) 
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for some positive constants vo and v\ provided that p and T are bigger than some fixed 
positive constants. Furthermore, K and © are compact operators on Ti. Therefore, us- 
ing the Fredholm alternative theorem we can conclude the existence of solutions to (4.17) 
provided that S n G V^H and T n + J(0 o ,/ n +i) G IC^H. These conditions can be veri- 
fied inductively, as in the usual one-component Boltzmann equation. We now write the 
conditions for n = which determine the macroscopic equations for p, u, T and ip: Since 
V[V v fo] and P[V v 0o] have no component along xo, it is easy to check that the condition 
VSo = can be written explicitly as 

d t p + V x • [pit] = 0, 

p[d t u + (u ■ V x )u] = -V X P + pK * p - <pK * <p, , (4.24) 
p[d t e + (u ■ V x )e] + PV X -u = 

where * denotes the usual convolution, P = pT is the equation of state for the pressure 
in the perfect gas and e = 3T/2 is its internal kinetic energy. On the other hand, since 
ICJ = and JCV V = 0, the condition ICT = becomes JCD t <j>o = which is explicitly 
written as 

d t <f + V x ■ [ipu] = 0, (4.25) 

Equations (4.24) and (4.25) represent the Euler equations for the binary mixture. They 
differ from the usual Euler equations by the presence of the equation (4.25) for ip and 
for the nonlinear self consistent force terms due to the long range Kac interaction. We 
will refer to them as the Vlasov-Euler equations (VE). Existence of solutions to the initial 
value problem for the system (4.24)-(4.25) requires some analysis but we do not discuss 
this. We simply assume that, for sufficiently smooth initial data a unique solution of the 
system exists and stays smooth up to some time t. 

Given such a solution, the functions fi and <p± can be found by solving (4.17) with 
n = 0. In consequence, fi is determined up to pi G VH and 0i up to q\ G KJH. The 
procedure can then continue by taking advantage of the arbitrariness of pi and q± to satisfy 
the conditions VSi = 0, KT\ = 0. In this way the functions f n and (fi n can be found for 
any n. Classical results by Grad [Gra] provide the smoothness and decay properties we 
use below. 
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Now we go back to the truncated expansions (4.7). Once <p n and n are computed for 
n = 0, . . . , K, we can look for the equations for the remainders Rf and R^. A straightfor- 
ward calculation shows that, in order that / and cj> satisfy (3.1), Rf and have to solve 
the equations 

D t R f + F ■ V v R f + W ■ V V R^ =e- 1 CR f + £ {1) R f + e m -\j{R f , R f ) + A f ], 

DtR<f> + F ■ V v R<p + W ■ V v R f =e- 1 TR 4> + e~ l QR f + (4.26) 



where 



and 



K 



C^g = J2 £h ~ 1 Wh,9) + J(9,f h )], 
h=i 

K K 

Qg = jC£e n <t>n,g), rMg = Jig^^fh), 

n=0 h=l 

l/=^" 2m+1 ( £ e h+h '- K -'Q(f h J h/ ) 

(h,h'): h,h'>l 
h+h'>K+l 

e h+h '- K [F h ■ V v f h . + W h • V v <f> h >] - A/k) 

(h,h'): h,h'>0 

h + h'>K 

(h,h'): h,h'>l 
h+h'>K+l 

£ e h+h '- K [W h ■ V v f h > + F h ■ V v <j> h >] - D^k) . 

(h,h'): h,h'>0 
h + h' >K 



K 



F=J2^ F n + e m K®R f , 

n=0 
K 

W= £e n W n + e m K© R+. 



(4.27) 



(4.28) 



(4.29) 



n=0 

The expressions of Af and show that it is convenient to choose K > 2m — 1 in order 
to get them bounded as e — > 0. 
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The construction of the solution Rf, of (4.26) is obtained using a fixed point argu- 
ment to handle the nonlinear terms in the equations. 

In Appendix A we shall sketch the proof of the Proposition below, which extends the 
similar result proved for the one-component Boltzmann gas without long-range interactions 
in [Ca80] and [Lac]. We assume for simplicity periodic boundary conditions, namely Q is 
the 3-dimensional torus of unit side. The potential of the long range force is assumed C°°, 
non negative and of compact support. Of course such assumptions could be relaxed, but 
we will not try to examine the most general setup. Moreover we use the norm 



sup 

veR 3 



e - 2 (l + \v\ 2 y/ 2 \f(-,v)\ J (4.30) 



and |/| s is the Sobolev norm of order s. 

We will refer below to sufficiently smooth solutions to the VE equations meaning solu- 
tions which are in H s for some sufficiently large s and such that the inequalities 

Tq<T<T u Po<p±(p<Pi 

are verified for suitable positive constants T , Ti, p and p\. 

Proposition 4.1. Suppose that (p,u,T,<p) is a solution to the Vlasov-Euler equations 
(4-24), (4-25) sufficiently smooth in the time interval [0,t\. Then there are positive con- 
stants £o an d C such that, for e < Sq a unique classical solution to the system (4-26) with 
m > 4 exists and satisfies the bounds 

sup \\R f (- ,t)\\ aAs < Ce sup [P/(-,*)|U,*, 8 + ||^(-,£)IUa s ], 
te[o,t\ te[o,£] 

(4.31) 

SUp \\R<t>(- ,t)\\ aAs < Ce SUp [P/(-,*)||a,£ ia +p0(-,*)||a,^J, 

te[o,t\ te[o,t] 

- - def 

for any positive a < T/2, T = sup xen te[0 ^ T(x, t), £ > 3, s > 2. 

Corollary 4.2. Under the assumptions of Proposition 3.1, for e < Eq there is a smooth 
solution (/ t e ,(/)f) to the rescaled Vlasov-Boltzmann equations (3.1) and moreover, denoting 
by M t the Maxwellian with parameters evolving according to the Euler equations, it satisfies: 

sup [\\f £ t - M t \\ aAs + M - ^M t \\ aAs ] < Ce. 
o<t<t Pt 
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5. Navier Stokes correction. 



The Navier-Stokes corrections to the hydro dynamical equations on the Euler scale are 
usually obtained by means of a suitable resummation of the Hilbert series expansion called 
the Chapman- Enskog expansion. For our purposes it is convenient to look at a modified 
version of the expansion proposed by Caflisch [Ca87]. 

We use the notation: for n > 0, f n = Vf n , fn = V 1 - f n , 4> n = JC(p n , <p n = K,^4> n , 
F„ = K /„, W n = K <j) n . The terms in the expansions are given as follows: we set 
M s :=M(l,u,T); 

fo = pM s , (f> = <pM s , / 1 = 0, 1 = O (5.1) 



Lh = V ± D t fo + F ■ V v f + W ■ 



(5.2) 



I>! = -J(<hJl) + lC- 



D t <l> - e^V[D t h] + Wo ■ V v f + F ■ V„0 O 
P 



(5.3) 



D t (fo + eh) + F • V w (/ + eh) + W • V„(0 O + efa) 



(5.4) 



K. 



D t (4> + efa) + Wo ■ V v {fo + eh) + F ■ V„(^ + efa) 



(5.5) 



£f 2 = -2Q(h,fi) 



+ V 



_L 



D t (h + ef 2 ) + [F ■ V v h + F 1 ■ V v f + W • V w 0i + W 1 • V„0 O ] 



(5.6) 



T4> 2 = -<pJ(M a , f 2 ) - J(0i, /i) + K, 1 - \D t 4> n + ^V[D t h] + sD t 4> 2 
+ [Wo ■ V w /i + V^i • V„/ + ' V v 0i + F l ■ VM] 



(5.7) 



Dth+ Yl [Fh-Vvfh'+Wh-Vvfa,] 



{h,h'): h>0,h'>0 

h+h'=2 



(5.8) 
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K. 



D t <h + i Fh ' V ^h> + Wh ■ V -^'] 



= 



(5.9) 



(h,h'): h>0,h'>0 

h + h'=2 



Before giving conditions defining the higher order terms of the expansion let us comment 
about previous conditions: Equations (5.1)-(5.9) have to be solved in the order they are 
written: (5.2) and (5.3) are used to determine fi and then 0i in terms of the hydrody- 
namical parameters (p, u,T, p) and their derivatives; as a consequence (5.4) and (5.5) only 
involve (p, it, T, tp) and represent the hydrodynamical equations we are looking for. Be- 
cause of (5.1) fi and <j>i are completely determined. Then (5.6) and (5.7) can be solved 
to find fi and 02, depending only on /o, /i, 0o, 4>i and on the hydrodynamical parts of fi 
and 4>2- Finally, (5.8) and (5.9) are linear equations in the hydrodynamical part of fi and 
<pi which can be used to determine them. 

We notice that the term proportional to e in (5.3) has been included to avoid third 
order derivatives in the hydrodynamical equations. This is usually done in the standard 
Chapman-Enskog expansion by expanding the time derivatives in powers of e. 

For n > 2 we set: 



k,j>i 

k+j=n+l 



+ v- 



D t (f n + ef n+1 )+ Yl [F h -V v f h/ + W h -V v (f) h ,] 



(h,h'): h,h'>0 

h + h'=n 



(5.10) 



T(f> n+1 = -(pJ(M s ,f n+1 ) - Y J (<PkJj) 



k,j>i 

k+j=n+l 



+ jc- 



D t (f> n + eD t j) n+1 + [F h -V v( p h , +W h -V v f h ,} 



(h,h'y. h,h'>o 

h + h'=n 



(5.11) 



(h,h'): h>0,h'>0 

h+h'=n+l 



= 



(5.12) 



K 



D t </> n+1 + i F h- V «<^' + W h ■ Vvfh'] 



{h,h'): h>0,h'>0 

h+h' =n+l 



(5.13) 
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The procedure used for and 4>2 can be repeated to get f n and (p n for any n > 2. 

As for the Euler limit discussed in the previous section, instead of looking for the 
convergence of the expansion we consider its truncation (4.7) where the functions f n and 
4> n are computed according to the procedure just explained, but setting f n = 0, n = 
for n > K + 1. The remainders Rf and have to be solutions of the equations 
D t R f + F ■ V v R f + W ■ V V R^ =e- 1 CR f + £ {1) R f + e m ~ x [j(R f , R f ) + Af], 
DtR^ + W -V v R f +F -V V R^ =e- 1 [TR^ + ®R f ] + T^R 4> (5.14) 

+ e m - 1 [j(R^R f ) + A (i> )] 

where F and W are given by (4.29), and 6 are given by (4.27), while the 

expressions of Af and A^ are slightly different from (4.28), and are: 

A f = e K ~^{ £ h+h '~ K ~ 1 Q(fhJh>)- 

(h,h'): h,h'>l 
h+h'>K+l 

e h+h '- K [F h -V v f h , + W h -V v (/) h ,} (5.15) 

(h,h'): h,h'>0 

h+h'>K 

- Dj K - V[F ■ V v f K + W • V„0k]), 

{h,h'): h,h'>l 
h+h'>K+l 

e^'-KlWh-VJv+Fh-V^h,] (5.16) 

(h,h'): h,h'>0 

h + h' >K 



- D t (px - K[W Q ■ V v f K + F ■ W v( j> K }). 



We now exploit the equations (5.2)-(5.5) in order to obtain the Navier-Stokes equations 
for the binary mixture with long range forces. 

Note that in (5.2) the force terms do not contribute because V ± 'V V M S = 0. Therefore 
fi has the same expression as for the one component gas without self-interaction: recall 
that 

3 3 
i,j = l i=l 



V ± (D t f ) = M 



(5.17) 
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with 

and v = v — u. 
Therefore 



. 1 _ v . . „ / v 5 _\ £j 

Aij = -(UiUj - y = I — - - T\ — 



3 3 
i,j = l i=l 



(5.18) 



with -01 and ^2 non negative smooth functions of \v\. Moreover, V[W v fi] = because fi 
is orthogonal to the invariants. On the other hand (M s x a , V v (pi) = for a = 0, . . . , 3, 
because 0i is orthogonal to the constants. Furthermore 



J dv[F -V v {fo) + W Q -V v M = 0, 



y efoCfFo • V w (/ ) + • V^o = -F oP - W 
Hence the mass equation is 

d t p + V-{pu) = 0, 



and the momentum equation is 



pD%u + VP = -eV • a + pF + ¥>W 



with = d t + u ■ V and 



(5.19) 



(5.20) 



(5.21) 



(5.22) 



a 



1-3 



-v(djUi + diUj - -SijV ■ u) 



and 



In order to compute the equation for the energy and for the concentration we need the 
expression of 0i, which has to be computed using (5.3). This implies 



1 ir± 



0! = r _i /c 



D t (j> - e^V[D t h] + W ■ V v f + F ■ V^o - J(0o, /1) 



(5.23) 
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We are interested in computing the component of V v 4>\ along Mx4 i-e., after an integration 
by parts: 

(M X4 , V^i) = - J dvvfa. (5.24). 
Since J(0 O , A) = jAfi ~ r/i, we get 

r-V^c/i) = -r-^A - A = ^r-^CA/o) - A, 

after using (5.2) to get the second equality The second term does not contribute to (5.24) 
since A is orthogonal to Mv. Hence, by (5.18), 

J dvvY^J^h) = ^J dvvT- 1 [v-VT^(v 2 -5T)]; (5.25) 

the term proportional to Vw is odd in v and hence does not contribute to the previous 
expression. Now we compute 

K,^[D t {<pM s ) + pW ■ W V M S + <f>F ■ V V M S \ = 

<pD?(M a ) + ipv • VM S + M s v -Wip + pW ■ V V M S + <j>F ■ V V M S (5.26) 



We have 



<pD?(M a ) = ^<pD?T(v 2 - 3T] - V v M s <pD?u 



Since the first term is even in v, only the second term contributes to J dvvT^ in (5.25). 
By using the equation for the momentum we get 

D^u = - [-VP + pF + <pW] - - V • a (5.27) 
P P 

Moreover 

2 2 2 

T^VP -<pF- + pW + <pF] ■ V V M S = [^VP + W P ~^ ] • V V M S . 
P P P P 



Now 



-VP • ^M s + M s v ■ Vcp = {-tp^VT - ^Vp + Vcp] ■ M s v = 
p T T p 

[-^IvT + pV^-M^ 
T P 
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and 



-<p±v ■ VTM S + <pv ■ VM S = [~VT + - ^)}<PV • VTM S + |m s £ ® 5 • Vu 
= ^VT • M s 5[^ - A] + | Ms £ <g) 5 • Vu. 

(5.28) 

The last term in (5.28) does not contribute to (5.24). Since the first term in the r.h.s of 
(5.28) cancels out with the term in (5.25) (remember that M = pM s ) we have 



r.h.s (5.26) = Mv V- - ^-{p 2 - v? 2 ) - £^V • a 

p p z l p z 



(5.29) 



It is now easy to check that in the computation of the l.h.s. of (5.23) the last term of 
(5.29) is compensated by the term e^V[D t fi\. Collecting terms we get 



(M s X4,V„0i) = —D 



V- 



W 



p p*T 



(P 2 ~ <P 2 ) 



where 



D := I dvMviV^Vi 



(5.30) 



(5.31) 



The other terms appearing in the equation for the energy are computed in the standard 
way: 



with 



Moreover 



dvx4 



ip • V/i + F • V v /i + Wo • V„0i 







/ 



V- dvxdi = -V ■ [kVT] 7 



-I 



K = I dvip2Bf. 



J dv X 4F • V w /i = - j dvv- F /i = 

j dv X 4W • V^i = - J dvv- VMi = -W ■ DQ, 
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(5.32) 



with 



<p 1 



W(p 2 - ip 2 ) 



(5.33) 



P P 2 T 



Therefore, the equation for the energy is 




(5.34) 



Finally, to get the equation for the concentration we have to exploit the condition (5.5). 



The last term has already been computed in (5.30). Therefore, the equation for <p is 



where Q is given in (5.33). 

Recalling the definitions of F and Wo we finally get the Vlasov-Navier-Stokes equations 
(VNS) for a binary mixture given by (3.5). 

As for the Euler limit, the arguments in the Appendix A prove the following Proposition 
5.1 which holds under the same assumptions as before Proposition 3.1: periodic boundary 
conditions and smoothness of the long range potential. 

Proposition 5.1. Suppose that for e > small enough there is a solution (p s ,u £ ,T £ , ip e ) 
to the Vlasov-Navier-Stokes equations (3.5) sufficiently smooth in the time interval [0,t\ 
independent of e. Then there are positive constants Eq and C such that, for e < Eq an 
unique classical solution to the system (5.14) with m > 4 exists and satisfies the bounds 




(5.35) 




d t <*p + V • (<pu) = eV 



(DQ) 



(5.36) 




SUp || J R/(-,*)||a,£, a <C'e SUp [P/(-,t)||a,^+||^(-,t)||a,^], 



SUp \\R4,{- ,t)\\ aAs <Ce SUp [P/(-,*)||a,£,a+p0(-,*)||a,*,J, 

te[o,t] te[o,t] 



(5.37) 



for any a < T/2, T = sup xeCijte[0 ^ T £ (x, t), £ > 3, s > 2. 
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Corollary 5.2. Under the assumptions of Proposition 3.1, for e < Eq there is a smooth 
solution (ft,4>t) t° the rescaled Vlasov-Boltzmann equation (3.1) and moreover, denot- 
ing by M t the Maxwellian with parameters evolving according to the Vlasov-Navier-Stokes 
equations, it satisfies: 

sup [||/ t £ -M t - sh\\ aAs + M - ^M t - eMa,i,s < Ce2 Y 

0<t<t Pt 



6. Incompressible Navier-Stokes limit. 



In this Section we consider a different scaling limit such that one can get hydrodynamic 
equations with dissipative terms of order 1 instead of order s as in the previous Section. 
We choose also in this case 7 = e but we use the parabolic space time scaling 

t ^ & X ^ & X . 

After rescaling, eq.'s (2.11) become: 

dtf + e~ x v ■ W x f + e~ x F ■ V v f + e" 1 W ■ V v 4> = e~ 2 J(f, /), (6.1) 

d t (f> + e~ x v ■ V x cj> + e~ x F ■ V v cj> + e" 1 W ■ V v f = e~ 2 J(<p, /), 

with F and W given by (3.2). 

We shall solve (6.1) as in the Euler case in terms of a truncated Hilbert expansion of 
the form 

K 

f = Y,z n fn + e m R f , 

n=0 

6.2 

K 



n=0 
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with suitably chosen positive integers K and m, but in this case we choose in a different 
way the terms of order in the expansion. 

fo = M , 0o = 

where M is a Maxwellian with parameters p and T some fixed constants and u = 0. This 
implies the vanishing of the forces at the lowest order: Wo = F = 0. We remark that 
choosing the first order term in the expansion to be a global Maxwellian is essential to 
the incompressible limit set up. The choice 4>o = is made for simplifying the computa- 
tions. By plugging (6.2) in the rescaled Boltzmann equations (6.1), one easily obtains the 
conditions for the higher order terms in the expansions: 

Ch = 0, (6.3) 

r>i = o, 

for < n < K - 1: 

cf n+2 + Q(fh,M = 

(h,h'): h,h'>l 
h+h'=n+2 

[ F h ■ Vvfh> + W h ■ V v (f> h ,] + V ■ V x f n+1 + d t f n , 

(h,h'): h>0,h'>0 

h + h'=n + l 

T(j) n+ 2+ ^2 J((f>h,fh') 

(h,h'): h,h'>l 
h+h'=n+2 

= \- Fh ' Vu ^' + Wh ' V ^'] + V ' V x</>n+l + QtK- 

(h,h'): h,h'>0 

h+h' =n + l 

Moreover 

d t R f +e~ 1 [vRf + F- V v R f + W • V„ify] 

= e~ 2 CR f + e^C^Rf + £ {2) R f + £ m " 2 [J(i?/, R f ) + Af], 
OtRtf, + £ _1 [v ■ R<f> + F ■ VyRtj, + W ■ V v R f ] = e~ 2 TR<j, 

+ e^T^R^ + e^e^Rf + T^R^ + Q^R f + £ m " 2 [J(^, R f ) + 
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(6.4) 



(6.5) 



(6.6) 



(6.7) 



where 

K 

£ {1) g = J(h, g) + J{g, hi £ (2) g = £ e h - 2 [J(fa, g) + J(g, fa)], 

h=2 

K 

Q^g = J(<f> u g), Q^g = J(£ £ - 2 ^, g), 

n=2 

r {1) g = J(g,h), = fy- 2 ^,/,), 

h=2 

A f = e K ~^\ J2 £ h+h '- K - 2 Q(hJv) 

{h,h'): h,h'>l 
h+h'>K+2 

e^'-x-^Fh.Vvfv + Wh-VM'] 

(h,h'): h,h'>0 
h + h'>K + l 

- d t f K -i - v ■ V x f K - edf K ^j 

(h,h'): h,h'>l 
h+h'>K+2 

eM'-x-^Wh-Vvfo+Fh-Vvfo] 

{h,h'y. h,h'>o 

h + h' >K + 1 

- d t 4>K-i - v ■ V x 4>k - ed4>K^j , 
and F and W are given by 

K 

F=J2^Fn + e m K®R f , 

n=l 
K 

n=l 

with F n = K © /„, W n = K © cf> n . 

We now find the expressions of the first terms in the expansions /i, 0i, /2, 02- 
From (6.4) we get 

(j)! = <piM Q 
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and from (6.3) we have 



v v 2 — 3T 
h = M Q {v) ( Pl + Ul ■ ? + ^ 2 



where pi,ui,$i are to be determined as functions of By (6.5) with n = we obtain 

P(vV a! /i+F 1 .V t ,/o) = 0. 



But 



while 



P[u • VJj] = M 



^2 _ 3T 

1 + 2 y 2 jV.-ux + y. V x (T P i + p6i) 



V[F 1 .V v f Q ] = -M j i -F 1 . 



Hence we find the conditions 



Va, • ui = 0, V x 



Pi + 1 + J dyU(\x - y\)pi(y)] = 0. (6.11) 



On the other hand, 



V ± [F 1 -V v f }=0, 



so 



-l 



-£-ij(/ 1 ,/ 1 ) + / 2 , 



with / 2 e NulLC. 

Therefore / 2 has the usual expression 



i,j = l i 1 



(6.12) 



^From (6.6) with n = we get the expression for 2 : 



2 = r- 1 v ■ v x (t>! + w v v f - J{<t>u h) 



+ 02, 



with 02 e Nulir. 
Moreover, by (6.3) 



-(fiT ^/i = -yi/i. 
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Hence 

<p2 = (flfl + 

^From (6.5) with n = 1 we get 



r- 1 [vM ]+4> 2 . 



and 



d t ui + ui ■ V x ui = -V x p + F 2 + vA x u x + piFi + tp{W\ 



-[dtOt + u x ■ VM =u 1 -F 1 + kA x 6± 



Since F 2 = V X G, with G(x) = J dyU{\x — y\)p2(y), putting p = p — G, the previous 
equation reduces to 

d t u\ + ui ■ V x ui = -V x p + vA x m + piFi +ipiWx, 

which is the usual incompressible Navier-Stokes equation with the self-consistent force 

piF 1 + (p 1 W 1 = -piV x J dyU(\x - y\)pi{y) + <piV x J dyU(\x - y\)(fi(y). 

Finally from (6.4) with n = 1 we get the equation for <pi 

rl If 

d t <pi + ut ■ V x (fi = D -A x Lp 1 - j;A x J dyU(\x - y\)(fi{y)] 



with 



D = -J dvv-T-^vM). 
Summarizing, dropping the index 1, the set of equations for p, u, 6>, 0, p is: 



(6.13) 



d t u + u- V x u = -V x p + vA x u + pF + (pW, 
-[d t 9 + u ■ VJ] = u-F + kA x 6, 



d t (p + u ■ V x (f = D 



1 1 f 

-A x ip-^A X J dyU(\x -y\)(p(y)] 
F = -V x J dyU(\x - y\)p(y), W = V X J dyU(\x - y\)cp(y), 
p + 6 + J dyU(\x-y\)p(y)] =0, 



(6.14) 



V x • u = 0. 
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The equation for <p is linear unlike the one we get in the VNS equations, but there is still a 
non linear term in ip in the momentum equation. The equation for 9 which corresponds to 
the deviation in the temperature decouples from the rest. In fact, if we consider a solution 
to the previous equation with an initial datum p = const, 9 = const such conditions persist 
in time and u and <p have to solve the simplified set of equations 



d t u + u • V ' x u = —V x p + uA x u + LpW, 



d t Lp + u- V x (f = D ^_A x Lp - i A x 
W = V X [ dyU(\x-y\)<p(y), 



dyU{\x-y\)(p{y)] 



(6.15) 



V T -u = 0. 



In Appendix B there is the proof of the following proposition: 
Proposition 6.1. Suppose that for e > small enough there is a solution (p, -u, T, tp) to 

the incompressible Navier-Stokes equations (6. 14) sufficiently smooth in the time interval 
[0,t\ independent of e. Then there are positive constants eo and C such that, for e < Eq a 
unique classical solution to the system (6.7) exists and satisfies the bounds 

sup \\R f (- ,t)\\ aAs < Ce sup [P/(-,*)||a,* l8 + ||^(-,£)||aA S ], 
te[o,F] te[o,t\ 

(6.16) 

sup \\R<t>(- ,t)\\ aAs < Ce sup [||A/(-,t)|| a> £ iS +||^(-,t)|| a) 4j, 
te[o,£] te[o,t] 

for any a < T/2, i > 3, s < m. 

Corollary 6.2. Under the assumptions of Proposition 3.1, there is for e < e$ a smooth 
solution (ff,(/)f) to the rescaled Vlasov-Boltzmann equation (6.1) and moreover, denoting 
by M t the Maxwellian with parameters evolving according to the incompressible Navier- 
Stokes equations (6. 14), it satisfies: 

sup | \f t £ - M - eh \\ aAs + | |0f - M - £0i | \ aAa < Ce 2 ] 
o<t<t 
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Appendix A 

We present a sketch of the proof of Proposition 3.1 in the case of hard spheres, where 
v(v) w \v\ for large v : s. The extension to more general cross sections is possible along the 
lines proposed in [DE], but we do not discuss it. The smoothness and decay properties 
of the expansion terms are obtained by now standard methods ([Ca80], [DEL], [ELM94], 
[ELM95], [ELM98], [ELM99]) which allow to prove the following 

Theorem A.l 

Given e > 0, assume that there exists a sufficiently smooth solution of the Vlasov-Euler 
equations (4-24) o,nd (4-25) in the time interval [0,t\. Then for any j > 3, s < so and 
< a < T* = sup^ ^gQxjQ ^ T(x,t) there is a constant c > such that the terms in the 
expansion fi, fa, i = 1, K, with K = 2m, solutions of the equations (5.8)-(5.7) 

\\fi\\a,j,8 < C, \\fa\\ a ,j,s<C (A.l) 
\\dvfi\\a,j,2 < C, \\d v fa\\ a ,j,2 < C (A. 2) 



We need bounds on the remainders Rf and R^. 
We set 

R r = RW = R f - R^ R b = R^ = R f + R^. 



(A3) 



The reason is that terms like W ■ V„-R^ and F ■ V v Rf in (4.26) are not well suited 
when some force term is present to represent the solutions of the equations in terms of 
characteristics, which is essential in the method we are going to use. The equations for 
the new variables are: 



+ r 2 i?« + +6i#( 2 ) + e m ~ l [J(RW,RW) + J(RW, RW) + A™] , 



D t R {2) + F {2) ■ V V R {2) 



P^jCrW + p^QR^ + pVfRW\ + C 2 R™ 
+ T X R^ + 2J R« + e m ~ x [J(RW,RW) + J(R (2 \ RW) + A™] , 
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(A4) 



where M is the Maxwellian M with p = 1 and 



£/i = J(M,/i) + J(/i, M), 

Oh = J(M,h), Th = J(h,M) 



K 



K 



3=1 3=1 

A^ = A f + A^ A^> = Af — Atf, 

Following Caflisch [Ca80] we now decompose the remainders in low and high velocity- 
parts, by looking for solutions to equations (A. 4) in the form 

RW = y/pWMgW + R {2) = y/pWMg™ + v / M^/i (2) 

M* is a global Maxwellian with a temperature T*. We have 

Dtg (i) + F (D . Vg W = £ -i UQLgW + y/pO)y/pl*)Tg<n + p^Gg™ 



VpW(K*/i (1) + K* T h^) + ±=K* G hW 

VP (1) 



Dth (D + F (l) . Vh (l) = a[li (D + F (D . ^]^pWg^ + F (D . ^ h (D 

r n( 2 ) i 

L p^ 1 - 1 J 



+ e 



m—l 



i/*Q*MV 1) + ^ (1) ^P (1 V 1) + ^ (1) ) 



The equation for g( 2 \ is obtained by the exchange 1 — > 2. 
where 

1, \v\ < 7 

X(v) = { X=l-X 
0, otherwise 
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(A5) 



(A.6) 



fj, 



(i) 



= ^A(logM), = \ P {i) V v \ogM,i= 1,2 



= -V^logM*, (j 



M 



(A7) 



If = -^—cVtif) = (-u + K)f, L*f = -L=CVJFf) = -u* + Kf 



M 



VM* 



ff = -^Q^/Mf, Qf = -Lfv^/ = (-vg + K G )f, , K* T 
VM 



1 



M 



G*f = —^T^fWf = (-v% + K G )f, df = - 1 L=T i VM*f, i = l,2 
LJ = -±=C t 4W*f, i=l,2, = ^=6,^/, z = l,2 



(A8) 



i/*Q*(/,*) = -^=Q(yfM^f,VM^e), v*J*(fJ) = -^=J(VM^f,VM^£) (A9) 



Linear problem 

We solve first the linear problem for gi,hi,i = 1,2, assuming that F w are given functions 
such that 

||F (i) ||oo + ||V :c F«|| 00 <« F . (A10) 

We consider the linear system 



,(2) 



1 VP {1) 



Dth W + F (i) . Vh (i) = a[fi (i) + F (i) . ^(i)] v^O)^(i) + F« • ^h™ 

,(2) 



+ £-V (1) 



- v + x(K*h^ + K* T h^) + ^{-va + K* G h^) 



+ £ m-l D (l) 

(All) 
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and the equation for g^ 2 \h^ obtained by the exchange 1 — > 2. Here has to be 
considered as a given force. 

We use the integral form of (A.ll) [ELM98]: 

9 {i) {t,x,v) = f dsH^(s^f_ t {x,v))e^{ -J*da'^(^)_ t (x,v))}. (A12) 



where i>W = p( l )i/ + p(^i/ Gj (p^'{x, v) the characteristics of the equation 

d t f + vV x f + F^-V v f = 



(A13) 



and 



+ 



e-V- 1 



,(2) 



VP (1) 



(A14) 



and if^ 2 ) is given by the same expression after the exchange 1 — > 2. 

v) = J* dsH'® (s, <p{i) s _ t {x, v)) exp { - J* ds'h® {<p(i) s ,_ t {x, v))}. (A15) 



with 



i>(0 = p(0 _ £ ^ . ir(0 



(which is positive for £ sufficiently small, depending on op, since grow linearly at high 
velocities) and 



+ L^iay^g^ + h™) + G^(ay/^)g^ + h™) + T^(a^g^ + h^) 
+ + F<V (1) ] y/^g^ + e m ~ l D^ 

(A16) 

We do not write explicitly the equations for g( 2 \ in integral form. In the following 
we use the compact notation: g = {g^\ g^} and h = {h^\ h^}. Below we use the 
notation || • ||^ s = || • \\qj jS and || • \\i = \\ ■ \\o.e,o- Generalizing the method by Caflisch 
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of g^ 1 ' , in terms of the L2 norm 



(A.17) 



[Ca80]to our case, we get bounds for the norms || 
of g^ 1 ' in the form 

||/i||r <£(l + a F )|M|L 2 +£ m |D|r-l 
Ibllr < ||<7||L 2 +£ m+1 |^|r, 

provided that £,£0 for some suitable £0 positive and finite for any finite To conclude 
the argument we need a bound for ||<7||l 2 = Yld=\ \ \9^\\l 2 m terms of the L2 norm of D. 
This last step is not standard so that we give a sketch of the proof. 

To estimate ||<7||l 2 , we multiply the first equation in (A. 11) by g( % \ i,j = l,2,z 7^ j, 
respectively, integrate over x, v and finally sum over % = 1, 2 



1 d 

2dt 



l<7 (1) ll! 2 + lb (2) ll! 2 



= s 



((v^g {1 \M- 1/2 J{M, V / M V / pW^ (2) )>+ 



+ (Vp^g {1) , M- 1/2 J(VMy/pWg {1) , M))+ 
(y/p^9 {2) , M" 1/2 J{ ^/My^g {2 \ M)) 



+ 

(2) 



(A18) 

Here (f,g) denotes the L 2 (0 x IR 3 ) scalar product. The operator L is symmetric with 
respect to (-,g • }. The terms in the first square brackets are non positive by the non 
positivity of the operator L. It is easy to see, by using the symmetry properties of the 
Boltzmann kernel [CC], that also the contribution coming from the terms in the second 
square bracket are non positive. Therefore we have 

1 d 



2 dt 



|^||| 2 <Ce- 1 ||/ i ( 2 ))|| r -|b|| L2 + /5|b|| 2 L 2 2 , 



(A19) 



with r > 3. The final estimate is 



\g\\ L2 < Ce m - 1 \\v- 1 D\ 
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The same kind of arguments provides the bound for the derivatives of g with respect to 
x. Because of the force terms the argument differs from the one given in [Ca80] in the 
fact that we have to control at the same time the derivatives with respect to velocity and 
space. We sketch the proof for the derivatives of g. Differentiating the first equation in 
(A. 11) we get two coupled equations for V v g and Vg 

d t Vg {l) + (v ■ V)V# W + VF« • V v g {l) + F« • V v (Vg {l) ) = VN&(g, h) 
dtVv9 (i) + „ . V(V^«) + V<7« + F« ■ V v (V v g^) = V v N^(g, h) 

where N(g, h) is the r.h.s. of (A. 11) for g. Proceeding as before in getting (A. 18) we obtain 
d 



dt 



(I|V^||l 2 + HV^IUJ < ca F (||V^|| L2 + HV^IUJ + ||VAT|| L2 + ||VJV|| L2 



The derivatives of N with respect to velocity can be estimated by the methods in 
[ELM94], where it is proven the identity 

9 )=w, |) + «c*§£) 

where stands for the partial derivative with respect to any of the components of v. The 
final result is 
Lemma A. 2 

There is an eo > finite for each finite cxf such that any solution to the linear problem 
(A. 11), with D and Fj given, satisfies for j > 3, s < 3 and any < e < Sq 



\V 1 D^\\j +2)S + \\U 1 U V U^ |j+2.,s 



d v D^\ 



(A20) 



H^II^ + H^Wll^-! < C{l + a F )e m - 8+1 \\\v- 1 D^\\^ 8 + \\v- 1 d v D 



)(0| 



\j,s 



(A.21) 



Non-linear problem and fixed point argument. 



The nonlinear equations (A. 5) are solved by a fixed point method. This method works 
if there is some small parameter in front of the non linear terms. In (A. 5) there are two 

37 



kinds of non linear terms: the usual Boltzmann non linear term, which is multiplied by a 
power of e and the Vlasov term involving the forces which gives rise to linear and quadratic 
terms in the remainders: in fact the forces are given by expressions of the type 

K 

= K © = k © [J2 £ n fi 2) ]) + £ m K © R {2) (A.22) 

n=0 

Hence the non linear term due to the force is small and we can apply the recursive argument. 

The Boltzmann terms are dealt with as in Caflisch [Ca80] . The control of the force term 
requires the boundedness of the gradient of the Kac potential. The result is 

Theorem A. 3 

There is an €q > such that the remainders satisfy for j > 3, z = s — (d — 1) with 
s < m and any < e < Eq 

\\9 {i) \\ j , z <ce m - 8+ \ \\h^\\ ] , z <ce m - s+ \ i = l,2 (A23) 



Proof. Let It£ be the solution of (A. 4) with force 



/V / (A.24) 



n=0 



:= + e m F ( k l) , j = i + 1 mod 2, 

and the collision integrals computed with replaced by R^-i', moreover, we put Rq = 
0,i = 1,2. Lemma A. 2 and an inductive argument assure that the sequences Rj£> are 
uniformly bounded for e sufficiently small. In fact, setting a = we have au = 

Q„(i) < a + C£ m ||.Rfc_i||i. By a standard argument ||i?fe_i||i is bounded by a constant 
\(ak-i) with A( • ) some monotone function. Hence, a = supctfc satisfies the inequality 
a < a + e m \(a). By setting e m \(2a) < a, we conclude a < 2a and hence the uniform 
boundedness of the sequences R^ for e sufficiently small. 

The differences 8Rf := R^ - Rf_ x can be decomposed again in high and low velocity 
parts 5gk,o~hk, with Sgk = (Sg^\ Sg^) and dhk = (Sh^ , Sh^ ) solve the equations 
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D t 5g k + (e m F k . 1 + F e ) ■ V v Sg k = N(5g k , 5h k ) 
D t 5h k + (e m F k . 1 + F £ ) ■ W v 5h k = N'(8g k , 5h k ) 
where N = (N (1 \ iV~( 2 )), N' = (N'^\ N'W) and N'W are the r.h.s. in (A.ll) with 



D« = s5F k ■ V v {g k -! + h k -!)+ 



[J(R k l v _ J ( R k-2i R kU)} + [ J ( R k-V R k-l) - J ( R k-2' R k-2)]- 

with SF^ := (F^ - F^\) = K © [R^ - R^_ 2 ). By Lemma A.2 the solutions satisfy 
(A.20) and (A.21), so that 



>(0 d(«) 



(0 o(j) 



>(0 d0') 



\\&9k\\j,s < Ce^-^Wu-^Fk) ■ V v (g k -i +hk- 1 )]\\ j+2 , 



\Sh k \\ jta < ^" s + 2 ||z>- 1 [5F fc • V w (^_i + 



(A25) 



(A26) 



We have 



\\SF k ■ V„((7 fc _i + h k -!)\\% < Csup(l + \v\ 2 y\V v (5g k + Sh k )\sup / dx\SF k \ 

x,v v J 

<C Jdx\Jdvj dyK(x-y)(R k . 1 -R k . 2 )^ 



<C dx 



sup(l + |w| 2 ) J / dyK(x-y)(5g k - 1 +5h k - 1 ) 



< C 



J dy\K{x-y)\ 2 ) 1 ' 2 snj>{l + \v\ 2 y{ J dy\5g k - X + 5/i fc _x| 2 ) 1/2 



<C[\\6g k ^\\% + \\8h k ^\\%] 

To get the second inequality we have used that the norms 1 1 V v Sg k \ \ j jS and ||V„5/ifc||j )S are 
finite by Lemma A.2, so that the supremum over x in the first row exists finite. The last 
inequality is a consequence of the fact that the Kac potential is bounded and that the 
space integration is on a torus. Finally, by (A. 25) and (A. 26) we have 

\\8gk\\j,o < ce m ~ S+1 \\5g k -i\\ JyS + ||^fc-l||j,o] 
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We remark that it is possible to prove that the norm || • \ \j jZ for the remainders g^ l \ hS 1 ^ 
are bounded with z = s — (d — 1) and s < m. □ 



Appendix B 



In this Appendix we show how to bound the remainders which are solutions of (6.7). 
The method we use is different from the one in Appendix A. In fact in this case the lowest 
order is a global Maxwellian and we do not need to introduce the decomposition into low 
and high velocity. Also in this case we need a Theorem on the regularity of the terms of 
the expansion analogous to Theorem A.l 

Theorem B.l 

Given e > 0, assume that there exists a sufficiently smooth solution of the incompressible 
Navier-Stokes equations (5.28) in (0,to]- Then there is a constant c > and s depending 
on the smoothness of the solution of INS such that the terms in the expansion fi,(f>i, i = 
2,...,K solutions of the equations (6.5) and (6.6) satisfy 

\\fi\\j,s < C, ||0i||j,2 < C 



\d v f, 



< c, 



\d v c 



< c 



(B.l) 
(B.2) 



for any j . 

We write (6.7) for the variables 1 = 1,2 defined in (A.3) 

dtR W + e -i [ v . R W + F (i) . V V R^} = e~ 2 [CRM + TRW + ®R(V 

+ e- 1 \dRW + T 2 R^ + <S> lR W] + \c[RW + r' 2 R^ + Q[R( 
+ e m ~ 2 [J(rW,rV) + J(RW, R ( V) + A™] , 
d t R™ + e- 1 [v ■ R^ + F 2 ■ V V R^} = e~ 2 \CR^ + TR™ + OR^ 1 







(B.3) 



C 2 R {2) + TiR® + 2 R {1) + C' 2 R {2) + T[R^ + Q' 2 R^ 

+ e m ~ 2 [J(bP\rW) + J(R( 2 \ RW) + A^} , 

6, £, T are defined as 6, £, f in the list after (A. 4) after substituting the global Maxwellian 
M (p,f) to M. Finally, 

C i g = J(fi i \g) + J(gji% e i9 = J(<pP ,g), T i9 = J(g, /<*>), 
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K 



h=2 



K 



K 



n=2 



h=2 



The first step is to consider the linear problem associated to (B.3), namely to study (B.3) 
with the last terms DW := J(RM,RW) + J(RW,RW) + A&,i = l,2,i ^ j given and 
fixed, independent of R^\ R^ 2 \ Moreover remembering that the forces vanish to 
the lowest order in e, we assume that the norms of and their gradients are bounded 
by some constant eolf- The role of the constant op is similar to the one discussed in the 
previous appendix and we do not repeat the iterative argument in this case. We will just 
provide an estimate for the L 2 norm in (x,v) \\RW2 '■= H-R^lh + H-R^lh for the solution 
R = (R^\ R( 2 ^) of this problem, the rest of the argument being standard (see for example 
[ELM98] ) . We put = v / Mo* (i) so that 



£ -2 |" L ^(i) + + TW (2)1 + £ -i \ LiR (i) + Q 2 y(l) + Ti ^(2) 



+ 



(BA) 



where the relation between the old operators C, 0, T, Ci, ©j, C i, &[, Y\ and the new ones 
L, T, G, Li, Ti, Gi, L' T' G\ is of the form 



41 



It is easy to see that, setting ||^||| = ||* (1) ||| + ||* (2) ||| , 



2 eft' 



1*1 



i=l 



-2 



(fW.Mo^fMo, v 7 ^ * (2) )> + <* (1) , M - 1/2 J(v^o * (1) ,M ))+ 



<^ 2 \M~ 1/2 J(M , * (1) )) + <* (2) ,M " 1/2 J( v / Mo * (2) ,M )) 

+ e" 1 [<*W, Lx*W) + <^^, L 2 ^( 2 )) + (G 2 *W + T!*< 2 >))+ 

(^(C^+Ta^ 1 *))] + [(*« L^^ + ^^L^ 2 )) 
+ (*^, (G" 2 #« + T[^)) + (^ 2 \ (G'^W + To'tfW))" 
2 n(») 

(5.5) 

First of all, we observe that the terms in the second square bracket are non positive [CC]. To 
estimate the other terms we will use the strict negativity of the operator L (see (4.20)) on 
the space orthogonal to the collision invariants and the following estimate on the operator 
J(/, h) (see for example [GPS]): for any Maxwellian M and for any y e [—1, 1] 



/ 



dv- 



\J(VMf, VMh)\ 
vM 



< [ dvu\f\ 2 [ dvv\h\' 

JR 3 JR 3 



(5.6) 



This inequality and the bounds on the / n 's imply the following bounds: 



5>«L 4 *«)| <C || h 2 II * W h 2 \\ M~ 1/2 f^ 



1.7,8) 



i=l 



i=l 



(B.7) 



|(^*«L^«>| < C || || La || * || La || M~ 1/2 J2U hs (B.8) 

i=l n=2 

where g denotes the projection of a function g on the orthogonal to the invariant space of 
L, while g is the projection on the invariant space. Note that the presence of the product 
j| y/is^> || 2 1| * || 2 depends on the fact that Li and L\ are both orthogonal to the collision 
invariants. 
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Similar estimates hold for the terms involving the other operators. By using the bounds 
on the / n 's and n 's and after some algebra, the terms in the forth, fifth and sixth rows 
are bounded by 

C || v 7 ^* Ik || * ||l 2 
To bound the last term in the first square bracket of (B.5), we note that 



i=i 



2 

L 2 ' 



and we assume e so small that e 2 olf < 1/2. 

We integrate (B.5) in time between and to- With the notation ty t ( • ) = ^( • > t), we 
get 



to \\L 2 — 



1 



V^^t III, +C F || 



(B.9) 



+ Cie- 1 + 1) || \\l 2 II * t \\l 2 +e m - 2 || £>( • , *) ||| 2 



The first term in the second line is due to the bounds (B.7) and (B.8). Moreover 

1/2 1/2 

II -^o /i llj'.s an< ^ II /i are bounded by the regularity of the solutions of the 

macroscopic equations for < t < T and || M -1 ^ 2 J2n=2 fn ILs< C by Theorem B.l. 
Using the inequality 

— \x 2 + (ci£ _1 + c 2 )xy < (ci + c 2 e) 2 y 2 /4 

valid for any positive e, x, y. with a; = Hy^^l^, y = \\^\\ and suitable constants c\ and 
C2, we get (since ^(-, 0) = 0) 



*(-,*0) ||L< / 



*(•,*) \\l 2 + \\Mo 1/2 D(.,t) HI 



In conclusion, by the use of the Gronwall lemma, for e sufficiently small, we get: 

sup || *(•,£) || L2 < C(t ,a F ) sup || M- 1/2 D(-,t) \\ La (5.10) 

0<i<i t€(0,t o ] 
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The bounds for the Sobolev norm of higher order in x,v are obtained by studying the 
equations for the derivatives as explained in Appendix A. Finally, writing the equations 
for the remainders in the integral form and using the property of the linearized Boltzmann 
operators of improving the regularity in v, we get the analogous of Lemma A. 2 

Lemma B.2 

There is an s$ > such that any solution to the linear problem (B.4 ), with D and given, 
after choosing K = 2m satisfies for j > 3, s < sq and any < e < Eq 



\g® | \ jiS + \\d v g® | \j, s -i < e m - 2 C F \\ |£><*> | | J+2 , S + | {v^D® \ \ j+2j 



(B.ll) 



| \h® | \j, s + I \d v h® I | JiS _i < e m - 2 C F \\ |D W | \j, s + | Iv^dvD^ 



\j,s 



(B.12) 



The dependence on the force in the bounds (B.ll), (B.12) does not affect the argument 
given in Appendix A to solve the non-linear problem, because in the bounds for the 
constant Cp h will depend on the norm of Rk-i- 

Appendix C 

To show formally the convergence of the microscopic one particle distribution functions 
to the solution of the VBE in the Grad-Boltzmann limit, let us consider the hierarchy 
for the rescaled correlation functions rj r j b of j r particles of species r and jb of species 6, 
defined as 

v ■ (z r ■■■ z r -z b ■■■ z h • t) - Nrl Nbl f 

dz r jr+1 ■ ■ ■ dz r Nr dz b Jb+1 ■ ■ ■ dz^NiS-'ql, S^qiy^S^ql v\ ■ ■ • , S-'q^S-'r)) 

(C.l) 

where z a = (q a ,v a ) is the phase space point of a particle of species a and is the 
probability distribution solution of the Liouville equation 

N N a Nf} 

9r m fiN + V i ■ VtiVN ~ A t E Yl V €f ^(lef - ■ V V? VN = 0, 

i=l a^f3 i=\ j=\ 
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which is valid in Tn, i.e. where the hard spheres do not overlap. On the boundary of Tn 
we assume the boundary conditions 



ftN(£i,vi,---, 6v, v N ; r m ) = hn(£i,vi, ■ ■ ■ , v'i, •••,€j,v' j ,---, £, N , v N ; T m ) 



if 



16 - e 



where = — u[u • (vi — Vj)], v'a = Vj + uj[uj ■ (vi — Vj)] with u> the unit vector directed 
as & — £j. The above conditions merely state the conservation of the probability during 
an elastic collision. As pointed out before, contacts of more than two particles have null 
Lebesgue measure, so they do not affect above definition. 

The rescaled correlation functions satisfy a hierarchy of equations of the form 

3b 



drr jr , jb + J2 [< ' V q rr jr , h + J2 V <^(K -Qj\)- V v rr jrtjb ] + 

i=l j=l 

3b jr 

E h 6 • ^ q f3r,3b + * 3 E v ^(i«? - • V wJ r irfc 

i=l i=l 

i=l 

[ B si r 3r,3b + l + ^;? r ir+l,ib + V ? br 3 r + h3b] , 



(C.2) 



+ 



i=l 



where, with the notation z fc = (zi, . . . , Zk), we have 

(^ + i.J(C;4) = 

/ dv r jr+1 / dw(v; r+1 -<) ■u\rj r +i, jb (g r + i y-,z% ) - r Jr+ltJb (z r jr+1 ; z b Jb )], 
Jr 3 Js\ 1 J 

Sl = {ue R 3 |H = l,w («; p+1 - <) > 0}, 

(^ r + l) = { Z ll • • • ( Z i) ■> ■ ■ ■ ■> Z j r l ( Z j r + l) )> 



with 
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■3r + l 



( r y r r -r \ 

— yz^, . . . z i , . . . , Zj r , Zj r +i) 



r Jr , Jb+ i((z r Jr Y; (z b Jb+1 )') - r jrJb+1 (z r jr ;z b jb+1 )j , 



and, for any z, z*, the phase points z' , z^, z and z*are defined by 

q' = <JS q'* = Q + 8v, v' = v — (v — v*), = + cu(u; ■ (t> — v*), 

g = g, q* = q — 8uj, v = v, v* = v*. 

Moreover, 

(^;[^^+i)(C;4) = 

where 

= • • • (A r ) 5 • • • > Zj r ), (z b j b+1 ) = (z b , . . . , Zj b , (Zj b+1 ) ), 

~~Zr~ I y r -r r \ b I y b b y b \ 

-jr ~ • • - Z i ' • • • ' ^'J' ^j b + l ~ ^1' • • • ' ^/'b> Z jb + lr 

The collision terms and are defined in a similar way. Furthermore 

(V?'V jVj6+1 )(4 r ;4) = 

- / K+i [ < + i7 3 v gr ^(kr-^ + ii)v, r ^ )ib+1 (4 r ;4 +1 ) 

and a similar expression for Vj' . Taking formally the limit 5^0, the limiting correlations 
satisfy the following Vlasov-Boltzmann hierarchy: 

jr jb 



i=l 



i=l 



+ B b > r r 



3r 

Yl \ B V r 3r + l,3 b ^ °i r 3r,3b + l 

1=1 

3b 

[ B i' br 3r,3 b +l + B ? br 3r + l,3b + Vl ,b r jr + 1 , jb ] , 



+ 



(C.3) 



i=l 



where, 



(^' r ^ + i )jb )(4 r ;4; 



'R 3 



dv 3r+i I duj ( vr j r +i - O • <*> r^+i^U^+i)';^) ~ r ^+i,i b (^v+i;4J 



46 



and, for any z, z*, the phase points z 1 ', z'^ are defined by 



q' = q'* = 9) v' = v — uj(u ■ (v — u*), = + cu(u; • (i> — u*). 

Similar modifications provide the other terms of the Vlasov-Boltzmann hierarchy. 
It is easy to see that if the initial condition of molecular chaos 



Jr Jb 



r jr,jb 



UL;4; ) = n/ r «'°)n/ 6 (^' ) 



i=l k=l 



is satisfied, then the correlation functions stay factorized at positive times r and f r (q, v, r) 
and f b (q, v, r) are the solutions of the coupled Vlasov-Boltzmann equations 



d T f r (q, v,t) + v VJ r (q, v, r) + F r ■ V v f r (q, v, r) = J(f r , f + 
d T f b (q, v,r) + v V q f b (q, v, r) + F b ■ V v f b (q, v, r) = J(f b , f + f b ), 



(C.4) 



where 



and 



F r (q,r) = -V q [ dq'^iVU^dq-q'l) f dvf b (q',v,r), (C.5) 

F b (q,r) = -V q [ dq'^iVU^dq-q'l) [ dvr(q',v,r), (C.6) 
Jo, Jtr 3 

J(f,g)= f dv, f du(v-v*)-u{f(v')g(v'J-f(v)g(v*)]. (C.7) 
Jtr 3 Js\ 

Summarizing, we obtained formally the Vlasov-Boltzmann equations for a binary mix- 
ture, where the Boltzmann collision kernel terms are due to the short range interaction 
while the Vlasov self consistent force is due to the repulsive weak long range interaction. 
If, instead of the hard core interaction the short range force is given by a finite range 
potential, we would get formally the same equations but with a different cross section. 

We want to stress that an important step is missing in order to make the above derivation 
rigorous. The first rigorous result on the derivation of the Boltzmann equation has been 
given by Lanford [Lan] where the convergence of the correlation functions is proven in 

47 



Loo-norms. On the other hand, the derivation of the Vlasov equation is based on the use 
of the variation norm and we have not been able to find a norm suited for both terms. 
The only related result, as far as we know, has been obtained in [GM] and is about a 
stochastic particle systems converging to a Vlasov-Boltzmann equation with a modified 
Boltzmann kernel (Povzner). The proof is based on martingale methods. In the linear 
case of a Lorentz gas with a Kac potential term it is possible to prove the convergence 
to a Boltzmann equation with a linear collision term and a non-linear self-consistent force 
term [MR]. 
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